^f- ■ 

O ' 

o : 

<N . 

> : 

From Diffusion to Anomalous Diffusion: A 
Century after Einstein's Brownian Motion 

5 1 

g : I.M. Sokolov 

Institut fiir Physik, Humboldt-Universitat zu Berlin, 
+^ ■ Newtonstr. 15, D-12489 Berlin, Germany 

and 
J. Klafter 

£h '■ School of Chemistry, Sackler Faculty of Exact Sciences, 

Tel Aviv University, Tel Aviv 69978, Israel 

February 2, 2008 

> ■ 

£NJ ■ Abstract 
CO 

C Einstein's explanation of Brownian motion provided one of the 

cornerstones which underlie the modern approaches to stochastic pro- 
^-j- ■ cesses. His approach is based on a random walk picture and is valid for 

, Markovian processes lacking long-term memory. The coarse-grained 

behavior of such processes is described by the diffusion equation. How- 
ever, many natural processes do not possess the Markovian property 
i . and exhibit to anomalous diffusion. We consider here the case of 

subdiffusive processes, which are semi-Markovian and correspond to 
q ' continuous-time random walks in which the waiting time for a step 

O ■ is given by a probability distribution with a diverging mean value. 

\ Such a process can be considered as a process subordinated to nor- 

JJp | mal diffusion under operational time which depends on this patholog- 

j_j ■ ical waiting-time distribution. We derive two different but equivalent 

forms of kinetic equations, which reduce to know fractional diffusion 
or Fokker-Planck equations for waiting-time distributions following a 
power-law. For waiting time distributions which are not pure power 
laws one or the other form of the kinetic equation is advantageous, 
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depending on whether the process slows down or accelerates in the 
course of time. 

The description of Brownian motion in Einstein's 1905 work re- 
lies on the assumption of the existence of a time-interval r, such 
that the particle's motion during different r-intervals is indepen- 
dent. The coarse-grained version of this motion leads then to the 
known diffusion equation. However, in many cases this assump- 
tion is violated. An example is dispersive transport in disordered 
systems which stems from a broad distribution of waiting times 
which may have a diverging mean. This ill-defined mean wait- 
ing time results in subdiffusion. In this contribution we derive, 
within a unified scheme, two equivalent forms of kinetic equations 
for subdiffusive behavior. For power-law waiting-time distribu- 
tions, the equations reduce to the "normal" form of a fractional 
Fokker-Planck equation with a fractional derivative replacing the 
first-order time-derivative, or to a "modified" form. For waiting 
time distributions which are not pure power laws one or the other 
form of the kinetic equation are shown to be advantageous, de- 
pending on whether the process slows down or accelerates in the 
course of time. 

In his epochal 1905 work, see [1], Einstein showed how the postulates 
of the kinetic theory of heat led to the conclusion that small but macro- 
scopic particles suspended in a fluid must perform an unceasing motion. He 
obtained the laws governing this motion which he identified as Brownian 
motion. The motion in the absence of an external force was shown to be 
diffusive. The discussion of the stationary state in the gravitational field led 
to quantitative relations from which, for example, Avogadro's number could 
be obtained with a precision which at that time couldn't be reached by other 
methods. The results of this and subsequent works [1] made up Einstein's 
PhD thesis defended in Zurich in January 1906. The approach proposed by 
Einstein included some of the modern approaches to stochastic processes, and 
stimulated and encouraged other scientists such as Smoluchowski, Langevin 
and Planck who had their own interests in molecular motion. Thus, Smolu- 
chowski in 1906 formulated a clearer probabilistic approach to the problem 
of diffusion based on what we call now random walks, which, according to 
his own words, was more straightforward and therefore simpler than Ein- 
stein's. Langevin in 1908 introduced an approach based on what nowadays 
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is called stochastic differential (or Langevin) equation, which was claimed to 
be "infinitely simpler". We should also mention that the widely used term 
" random walk" stemmed from a question put forward to readers of " Nature" 
by Carl Pearson in 1905 motivated by a biological problem. 

Einstein's description of the motion of suspended particles was based on 
three postulates: 

• The motion of different particles is independent (which is always true 
if their concentration is low enough), so that the problem is essentially 
a one-particle one. 

• There exists a time interval r such that the displacements of the same 
particle during different r-intervals can be considered as independent 
(fast decorrelation). 

• There exists a mean squared displacement A 2 of the particle during 
such a r-interval . 

The second postulate guarantees the Markovian nature of displacement 
at times larger than r, and the third one leads to the convergence of the 
corresponding process to a normal diffusion (and, say, not to a Levy flight). 

The overall picture put forward by Einstein looks as follows: The state 
of the system (the coordinates of the Brownian particle) is sampled at time 
intervals of increment r and the displacement during each interval r is chosen 
according to some probability density function (pdf) p(s). This picture is 
essentially a random walk picture since considering the system only at times 
which are the multiples of r leads to the notion of steps: the pdf p(s) plays 
the role of the distribution of step lengths, and the waiting time until the 
next step is exactly r. 

The formulation in terms of discrete jumps or steps is adequate in solids, 
where one of the mechanisms of conduction is attributed to thermal activa- 
tion out of bound impurity states. However, here the idea of a fixed waiting 
time r between the jumps has to be abandoned, and the waiting-time dis- 
tribution ip(r) of the times between subsequent steps has to be explicitly 
taken into account. This defines the continuous time random walk process 
(CTRW) which was first introduced into mathematical physics by Montroll 
and Weiss [2], and applied to semiconductors in a seminal work by Scher and 
Montroll [3]. The decoupled CTRW scheme assumes that the step lengths 
and the waiting times are independent random variables, chosen from cor- 
responding distributions with probability densities p(s) and ip{r). A simple 
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argumentation based on an exponential form of the tail of density of im- 
purity states and on the Arrhenius law leads to the power-law form of the 
waiting-time distribution [4], 



tp(r) oc r 



-l-a 
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where the exponent a is proportional to the temperature. Of particular 
interest is the case of < a < 1 for which all moments of ip{r) diverge. This 
process has no characteristic time scale, a strong violation of the original 
1905 picture, and results in subdiffusion [5-8]. The Markovian nature of 
the process is violated as well, and replaced by a semi-Markovian character 
[9]. In what follows we use the scalar notation; however, no problems arise 
when taking x to denote the components of the particle's position in three- 
dimensional space. 

We note that the particle's displacement under a random walk process is 
given by the stochastic differential equation: 



where a* gives the length and direction of the corresponding step taking 
place at the time instant ti. Coarse-graining this picture over some typical 
time interval At leads to a picture which corresponds to the overdamped 
Langevin dynamics: if both the mean time between steps (r) = — 1{) 
and the mean squared step length (af) exist, then averaging over the interval 
At >> (r) leads to appearance of a Gaussian noise with intensity (aj) jr. 
In this coarse-grained sense one can also put down 



The noise is uncorrelated on intervals At ^> (r). If the mean waiting time 
diverges, < a < 1 in Eq.((TJ), no convergence to a Gaussian takes place for 
whatever long t due to the existence of indefinitely long step-free intervals. 
The Langevin scheme for CTRWs with power-law waiting time distribution 
assumes thus the introduction of correlated non-Gaussian noise [10] and loses 
its appealing properties as intuitive description instrument. 

The original approach to CTRW based on the Laplace-Fourier-represen- 
tation of the process can be generalized to models showing correlations (like 
Levy walks [11]) but, due to the crucial assumption of spatial homogeneity, 
can hardly be used for the description of systems in external fields unless 
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these are homogeneous. On the other hand, the uncorrelated CTRWs can 
be considered easily also for spatially inhomogeneous situations, as long as 
this inhomogeneity concerns only step lengths, and not waiting times. This 
is so, since this latter case corresponds exactly to the case of subordinated 
random processes: if almost all t/>(t) (except at most the ones corresponding 
to boundary sites) are the same, the dynamics of the system can be considered 
as a process subordinated to simple random walks. To elucidate the situation 
let us start from the system without physical boundaries. 

Considering the situation not as a function of time, but rather as a func- 
tion of the number of steps, we can easily convince ourselves that the par- 
ticle's displacement as a function of the number of steps is a discrete time 
random walk, possibly with position-dependent step lengths. Let us con- 
sider the number of steps as the internal, operational time governing the 
system's evolution. In this operational time the evolution is given by the 
x-dependent probabilities for discrete jumps, or by probability densities for 
continuous step lengths. In what follows we discuss the continuous situation, 
and describe the displacements by their pdf's; the discrete case differs only 
terminologically. The role of the waiting times in this case reduces to the fact 
that the actual number of steps made up to the time instant t fluctuates, so 
that the operational time is a random function of the physical time t. The 
fact that this function is monotonously nondecaying with t and thus allows 
for causal ordering of the events (i.e. is really a time), is in general impor- 
tant, but plays no role in the following considerations. From the probabilistic 
point of view, the distributions for subordinated processes at some physical 
time t are mixtures of the corresponding distributions of the underlying pro- 
cess at different operational times. For example, for the pdfs of processes 
with continuous distributions of step lengths one has 

P(x,t) = J2w(x,n) Xn (t) (3) 

n 

where W(x, n) is a probability distribution to find a random walker at point 
x after n steps, and Xn{t) is the probability to make exactly n steps up 
to time t. In a classical Scher-Montroll CTRW for which the waiting time 
distribution follows Eq.((TJ) with < a < 1, Xn{t) corresponds to a random 
process in which n typically grows sublinearly in t, n(t) ~ t a . One can 
say that the operational time is always in delay compared with the physical 
one, see Fig.l. Thus, the overall process is subdiffusive. We note that some 
superdiffusive processes can also be considered as subordinated to random 
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walk; in this case the typical value of n grows faster than linear in t [12]. In 
the present work we concentrate on the subdiffusive case. 




Figure 1: The operational time of a CTRW process with a power- law waiting 
time distribution between steps. The steps follow very irregularily in physical 
time, and show step-free intervals on all time scales. On the average, the 
operational time is in delay compared with the physical one. 



In the continuous limit one can pass from a sum to an integral, changing 
the discrete variable n to a continuous variable r, In this case 

POO 

P(x,t)= W(x,r)T(r,t)dr. (4) 
Jo 

The function T(r, t) for a decoupled subdiffusive CTRW can be easily 
calculated, since the probabilities Xnif) are known [2,6]. Their Laplace trans- 
forms are given by 

1-jiCu) 

XoW = 

u 

and 

Xn{u) = Mu)r-\u)^—^- (n > 1). (5) 

u 

Here and in what follows the Laplace-transform of a function of time is 
denoted by tilde and ipi(u) is a Laplace transform of the forward waiting 
time density of the first step ipi(t), which might differ from ijj(t) for all 
subsequent steps. In what follows we take that the beginning of the time 
count coincides with the preparation of the system ("zeroth" step), so that 
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■01 (u) = ip(u). Moreover, this corresponds to the assumption that for t — 
one also has r = so that the initial conditions for probability densities 
P(x,t) and W(x,t) coincide. For any n we thus have [2,6] 



Xn{U) 



u 



u 



exp 



-n\ni/j(u) 



(6) 



Since ip{u) is a Laplace-transform of a pdf, it is continuous. In addition 
ip(0) = 1, which allows us to introduce the function <f>{u) = 1 — ip{u) which 
is small for u small (corresponding to large t). In the time domain <ft(t) is 
simply given by <fi(t) = 5(t) — ip(t). Expanding the logarithm we obtain 



f(n, u) 



XnW 



4>(u) 



U 



exp 



-n0{u) 



(7) 



Applying Laplace transform to both sides of Eq.(j3j), 



P(x, u) 



W{x,T)f(r,u)dT 







W(X,T 

o ' u 

(u) ^ 



u 



exp 



-T0[U 



U 



W(x, t) exp —T(f>(u) 



dr 
dr 



u 



(8) 



i.e. P(x,u) can be obtained from W(x,u) through a change of variable. 
The particular case of <p(u) ~ u a is considered in [13]. One can make a 
step further and obtain the equation governing the behavior of P(x, t), if the 
equation giving the "time" -evolution of W(x,t) is known. Let us assume 
that W(x, t) satisfies the following equation 



■%-W(x,t) = £ x W(x,t) 

CLT 



(9) 



where the linear operator C x acting on the x- variable of W might be discrete 
as in the master equation, or continuous as in the diffusion equation with 



dx 2 
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or in a Fokker-Planck equation 

d 2 d 

£x = £-FPE = K—r - fj,f(x) — , 

OX ox 

where K is the diffusion coefficient, fi is the mobility, and f(x) the external 
force. It is only important that it does not depend on any of the time variables 
of the problem, t or r (the last assumption can be relaxed, see [14]). The 
Laplace-transform of Eq. © reads 

sW(x, s) - W(x, 0) = C x W(x, s) (10) 

where s is the Laplace variable. Taking s = 4>{u) we get 

4>(u)W(x, 4>(u)) - W(x, 0) = C x W(x, 4>(u)), (11) 

and multiplying both sides by <p{u)/u 1 we obtain the equation for P(x,u) 

4>(u)P(x, u) - ^^W(x, 0) = C x P(x, u). (12) 
u 

From this equation different forms of the corresponding kinetic equation in 
the time domain can be obtained. These forms are essentially equivalent, 
however, depending on the analytic properties of the function <f)(u), one form 
or another can be preferable for practical applications. 

Returning to the time-domain we get from Eq. ( [T2*|) the following integro- 
differential equation: 

f (f>(t-t')P(x,t')dt' -W(x,0) [ <j)(t')dt' = £ x P(x,t), (13) 
Jo Jo 

where the function <f>(t) is the inverse Laplace transform of the function <f>(u). 
Inserting <ft(t) = 5(t) —i/j(t), performing partial integration in the convolution 
term, and using the fact that P(x, 0) = W(x, 0), one arrives at 

C(t - t')^;P(x, t')dt' = C x P(x, t). (14) 



^^The function ip(u) is a Laplace transform of a pdf and therefore is completely 
monotonous. This means that it is positive and monotonously non-growing. Moreover, 
ip(u) = 1 so that 4>{u) is monotonously nondecaying and (j>(0) = 0. For all distributions 
except for S(t) it is monotonously growing and thus possesses no zeros except for one for 
u = (which corresponds to t = oo). Therefore all divisions or multiplications discussed 
in the text are essentially harmless operations. 
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The integral kernel C(t) here is connected to the cumulative distribution 
function of waiting times, C(t) = 1 — F(t) = J°° ?p(t')dt' , i.e. it is the 
probability to make no step up to time t. This form of the kinetic equation 
will be termed here as a generalized Caputo form, for the reasons which will 
be clear below. 

The kinetic equation, Eq. (jl4|) . has a clear mathematical meaning, and 
nonnegative solutions (subordinated to those of the Markovian evolution 
equation, Eq.()10|)) for all nonnegative, monotonously decaying kernels C(t). 
Essentially one can even say that all equations of the form of Eq. ()14|) whose 
kernels are monotonously non-growing, nonnegative functions describe con- 
tinuous limits of decoupled CTRW. 

A simple example of the exponential waiting time pdf 

V>(t) = Aexp(-Ar) (15) 

leads to C(t) = exp(— At). For t ^> A -1 the kernel C(t) may be approximated 
by a 5-function, C{t) ~ A _1 5(t), and the integration over t' leads to 

^P(x,t) = \jC x P(x,t), (16) 

corresponding to a simple time-scale change by a factor of A. The process is 
back to Markovian. 

In case of power-law distributions ip(t), Eq.((TJ), the left hand side (lhs) 
of Eq. (jl4J) defines [16,17] an interesting mathematical object, namely a frac- 
tional Caputo derivative. The Caputo derivative is an operator of the form 

= = f(ibj jf <** " T) ^ /(T) - (17) 

Here o/f is a fractional integral operator defined through 

oI?g(t) = f ^ ) j\t-Ty- l g(T)dT, (18) 

< (3 < 1. The Laplace representation of this operator corresponds to mul- 
tiplication by . Thus, the Laplace representation of oD^ t is D^ t f(t) ==' 
u~ 1+a (uf(u) — f(0)), exactly the form in the lhs of Eq. (jl2J) for <f)(u) = u a . 
This equation corresponds to a "normal" form of the kinetic equation [15] 
with a fractional operator instead of first-order time-derivative on the lhs: 

D«P(x,t) = C x P(x,t). (19) 
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Integro-differential equations of the type of Eq. (|14|) do not have the form 
one typically uses. We now obtain a conjugated form, resembling the gen- 
eralized Fokker-Planck equation [18]. Multiplying both sides of Eq. (112)1 by 

u/(j)(u) we get 

uP{x,u) - W{x,0) = uJ—C x P(x,u), (20) 

4>(u) 

which is the Laplace-transform of a generalized kinetic equation in the form 

^P(x, *) = j t f ®(t ~ t')£ x P(x, t')dt'. (21) 

The lower integration limit follows from the fact that the Laplace representa- 
tion of the corresponding integro-differential operator coincides with u/<p(u) 
only if the integral term vanishes at the initial condition. The main differ- 
ence when compared to the typically used forms [19] is the additional time 
derivative in front of the integral. Note that the memory kernel $(t) of the 
integro-differential operator has a clear physical meaning. Let us expand 

J- = 1 = 1 + j(u) + 4> 2 (u) + ^(u) + ... (22) 
<p{u) 1 — ip(u) 

so that 

$(t) = S(t) +ip + ip* ip(t) +?p*ip* ip(t) + ... (23) 

where tp*...*ip(t) = J dt\... J dtn^ltijipfa)...^^ — ^) is the n-th convolution 
of the pdf ijj(t) with itself, which gives us the pdf that the n-th step of the 
CTRW takes place at time t. The sums of these probabilities gives us the 
density of steps: Q(t)dt is the mean number of steps taken in CTRW in the 
time interval between t and t + dt, which makes the overall picture given by 
the equation rather transparent. Thus, the integral kernels for the equations 
describing CTRWs have to be nonnegative. 

Contrary to Eq. (jl4j) . which is a valid kinetic equation for any nonnega- 
tive decaying kernel, the analytical properties of Eq. ()2T)) are more complex: 
some kernels do not correspond to any decoupled CTRW, others even do not 
warrant the non- negativity of solutions (i.e. are "dangerous", see Ref. [20]). 
Let us show that a necessary condition for Eq. ()2T|) to describe the decoupled 
CTRW is the divergence of the integral J °° Q(t)dt. To see this it is enough 
to note that for small u one has <p(u) — > 0, so that the Laplace transform 



10 



$(«) of the kernel <&(i) has to diverge for small u; otherwise the normaliza- 
tion of the waiting-time density is violated. This divergence means that the 
integral of the kernel over the time must also diverge. Equations with inte- 
grable kernels do not describe any decoupled CTRW-like process, although 
they might appear in a context of coupled motion [21]. As an example we 
see that the exponential kernel = exp(— t/t ) (for which Eq. (J21|) with 
a Fokker-Planck operator C x reduces to a kind of a telegrapher's equation) 
can never appear as a true kernel of the decoupled problem [20]. 

For a Poisson process characterized by the waiting-time pdf, Eq.(|15J) the 
density of steps tends to a constant $(t) — > A, so that applying the integro- 
differential operator is equivalent to a multiplication by a constant, which 
corresponds only to the change of the time-scale. 

For the power-law waiting-time density, Eq.((TJ), we have 4>(u) oc u a , so 
that oc u~ a so that oc t~ 1+a is a kernel of a fractional integral 

operator Q If defined by Eq.tflHjl. The integro-differential operator in Eq. (}2"Tj) 
corresponds to the operator of the Riemann-Liouville fractional derivative of 
the order 1 — a, see e.g. [7,8]. This is what we call the "modified" form of 
the fractional Fokker-Planck equation introducing an additional Riemann- 
Liouville fractional derivative on the right hand side (rhs) of the equation 
which now reads 

^P(x,t)= Q Dl- a C x P(x,t), (24) 
where the operator -Df is defined through: 

^^s^^fa^sjf^'-^w- (25) 

Pure fractional kinetic equations such as Eq. (jl9j) and Eq. (j24j) are not the 
only "reasonable" forms. In general both forms, the "normal" and the "mod- 
ified" one, are equivalent. However, for waiting-time distributions which 
asymptotically are not pure power-laws one or the other might be preferable 
for practical applications, depending on the properties of the waiting-time 
distribution. Special cases are kinetic equations with a distributed-order 
Caputo derivative [17,22] 

f daw x {a) D l - a P(x,t) = C x P(x,t), (26) 
Jo 
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and with a distributed-order Riemann-Liouville derivative [15] 



d_ 

di 



P(x,t)= [ daw 2 {a) Dl~ a C x P{x,t), (27) 
Jo 



where Wi(a) are some weight function. 

The "normal" form with distributed-order Caputo derivative is well fitted 
for describing processes getting more anomalous in the course of the time 
(retarding subdiffuion, as exemplified by crossover models of Refs. [17,22] 
and by Sinai-like diffusion processes [23]), while the "modified" form with a 
distributed-order Riemann-Liouville derivative describes the processes which 
in the course of the time get less anomalous (accelerating sub diffusion). Let 
us demonstrate these statements with two examples of different waiting- 
time distributions corresponding to crossover behaviors. Here it is easier to 
start from the asymptotic forms of cumulative distribution functions \l/(t) = 
l-C(*) = /o^(*')<ft'>ft> r *» 1: 

- 1 - £ - £ (28) 

and 

tt 2 (t) ~ 1 1 —— (29) 

y ' ct a + dtP y ' 

corresponding to an arithmetic and a harmonic means of the two power-laws 
with powers a and [3, a < (3. Both waiting-time distributions correspond 
to a crossover between two regimes each governed by a different power. In 
the case of short times are dominated by the larger power j3 and longer 
times are dominated by smaller power a. The crossover between the two 
regimes takes place at t ~ t c = (b/a) 1 ^ l3 ~ a \ In the case of ^ the behavior 
at short times is dominated by the smaller power a while the behavior at 
long times is dominated by the larger power (3. The crossover takes place at 
t c ~ (c/cf) 1 /^") . The long-time asymptotics of the corresponding waiting- 
time pdf's are 

7^1 + 7^1 (3°) 
which behaves as ^(0 — ^frr f° r t ^ *c an d as ^(t) — ^fr f° r * ^ *o an d 

cat"- 1 + d/?^- 1 
^ = (ct^dt^ (31) 



12 



o 



-10 



-14 




4 6 

log t 



10 



Figure 2: The probability density functions ipi(t), Eq. pOj) (full line), and 
ip2{t), Eq. (j3~T|) (dashed line). The parameters are a = 0.1, /3 = 0.9, a = <i = 

to- 3 , 6 = c = i - nr 3 . 



— for t < t c and as ^(A) 



- for t > t c . A 



which behaves as ^if) 
particular example of such pdf 's is given in Fig.2. 

Let us discuss first the types of the CTRW processes described by such 
waiting time distributions. To do this we discuss the mean number of steps 
(mean operational time) as a function of t. This quantity is of immediate 
importance since it is proportional to the mean squared displacement in a 
free diffusion or in a homogeneous external field. To obtain this value we 
note that from Eq.@ it follows that 



n[u) 



u 



i i 



n=0 



U 



u6(u) 



(32) 



which corresponds to the time-integral of the function $(«). Let us calculate 
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the functions <3>(w) for both our distributions. For the case of we get: 



^(u) = I - {Au a + Bu 13 ) (33) 
so that the function is 

= = i (34) 

U ; l-^i(u) Au a + Bu?> V ; 

and 

= Iau°IbuP (35) 
and use the Tauberian theorem, stating that the Laplace-transform of a 
function f(t) ~ t~^L(t), where L(t) is a slowly changing function of the 
time (for which lim^oo (L(kt)/L(t)) = 1 for any positive constant k) cor- 
responds in the Laplace representation for u — > to a function /(-u) ~ 
M /3-ir(l - f3)L(l/u). In our case 

= Au^ 1 1 + (B/A)uP- a = Au^i + ^B/A)^/^ - ^ 

is a product of a power-law and of a slowly changing function of the vari- 
able \ju and thus is a Laplace transform of a function which asymptotically 
behaves as 



Ar(a + 1) 1 + (5/A)^-/ 3 

= 1 1 (37) 

r(a + l) At~ a + Bt-P' y } 

In the case of ^2 the function $2(1*) = 1/02 (u) can be also obtained 
using the Tauberian Theorem. To do this we note that the function 2 (t) is 
a derivative of C^i), and therefore its Laplace-transform is 

<h(u) = uC 2 (u). (38) 

On the other hand, C^w) is the Laplace transform of the function 

° 2{t) = ^Tdt? = If (i + ( c /d)^) (39) 
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being a product of a power-law and a slowly changing function. Therefore 

c 2 (u) ~ r(i - ^crV-^i/u) = r(i - ^'V- ^J^^ . (40) 

The Laplace transforms, $2(w) and 77,2 (w) for this case read 

$ 2 (n) ~ i$ 2 («) = 1 [d^ 1 ^ + en" 1 -] (41) 

and 

fiaCt*) ~ ±$ 2 (u) = 1 [du- 1 ^ + ciT 1 -"] , (42) 

so that 

i r d 

n 2 (*) 



r(i-/3) |r(i + /3)" 1 r(i + «)" J ' (43) 

Since all Gamma-functions in Eqs. (|37|) and (J4H)) are of the order of unity, the 
typical crossover times between the short- and long-time asymptotic behavior 
for n(t) are the same as in the corresponding waiting times. 

Let us now compare the forms of the corresponding generalized diffusion 
and (Fokker- Planck) equations. To obtain these generalized equations we 
have to calculate the kernels in Eqs. (jl4j) and ()21jl. For the first function, the 
behavior of C, follows trivially: 

Ci(*) = £ + £- (44) 

Taking into account the definitions of the corresponding fractional operators, 
we see that the generalized diffusion equation has a simple form with two 
fractional Caputo derivatives on the lhs: 

(A, B« + B l dQ P{x, t) = C x P(x, t) (45) 

with Ai = T(l — a) a and Bi = T(l — f3)b. Quite opposite, the function 
^2^) does not correspond to any simple form of the equation. The kernel 
C 2 {t) does not possess a form resembling any fractional operator, except 
in its far asymptotics C(t) oc t _/3 . This confirms our statement that the 
" normal" form, reducing in our case to the form with two Caputo derivatives 
of different orders, is advantageous when describing the processes which get 
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more and more anomalous (in our case: slower and slower) in the course of 
time. 

Let us now turn to the "modified" form. Using the Laplace-transforms 
of the corresponding functions, Eqs.tjHU) an d W^ 1 we see that $i(w) does 
not correspond to a Laplace transform of any simple fractional operator. On 
the other hand, the Laplace transform of the function $ 2 corresponds to the 
Laplace-transform of the sum of two kernels of Riemann-Liouville fractional 
integrals of the orders of a and (3. Therefore the corresponding "modified" 
equation has two Riemann-Liouville fractional derivatives on the rhs 

-P(x, t) = (A 2 D}- a + B 2 D]-^ C x P(x, t) (46) 

with A 2 = c/r(l — f$) and B 2 = d/T(l — j3). In general, the "modified" form 
is advantageous when describing processes which get less anomalous (here 
less subdiffusive) in the course of time. 

In summary, we have considered a decoupled CTRW model, in which the 
directions and the lenghts of steps, and the waiting times at sites are mutu- 
ally independent random variables. Such a CTRW can be considered as a 
random process subordinated to simple random walks. Using the idea of tem- 
poral subordination, we derive two different, but equivalent, forms of kinetic 
equations describing their long-time behavior. The corresponding equations 
the integro-differential ones, and they reduce to the fractional equations with 
a Caputo- and a Riemenn-Liouville derivatives for waiting-time distributions 
characterized by the asymptotic power-law pdfs. We discuss the advantages 
of each of the representation by considering two different pdfs showing a 
crossover from one power-law behavior to the other. The corresponding ki- 
netic equations might then take a form of fractional equations with Caputo 
or with Riemann-Liouville derivatives of distributed order. 

The authors are thankful to A.V. Chechkin, R. Metzler and R. Gorenflo 
for valuable discussions. IMS gratefully acknowledges the support by the 
Fonds der Chemischen Industrie. 
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